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摘 要 : 为 了 解决 零 水 印 算 法 存在 抗 几 何 攻击 能 力 弱 的 问题 ， 提 出 了 一 种 抗 几 何 旋 转 攻击 的 零 水 印 算 法 。 首 先 根据 尺 
度 不 变 特征 变换 (SIFT) 旋 转 校正 后 图 像 像素 失真 情况 ， 在 中 心地 带 确定 像素 近似 无 损 的 安全 区 域 ; 其 次 将 该 区 域 进 

二 级 宛 余 离散 小 波 变 换 提 取 低 频 区 域 ， 对 该 低频 区 域 进行 分 块 并 提取 每 个 块 的 最 大 奇异 值 ， en 
构建 过 渡 算 阵 ; 然后 通过 比较 过 渡 憩 阵 的 每 个 元 素 值 与 其 均值 关系 构造 特征 矩阵 ; 最 后 将 加 密 后 的 水 印 图 像 与 特征 矩 
阵 共 同 构建 零 水 印 。 实 验 结果 表明 ， 与 仅 利 用 SIFT 校正 算法 相 比 ， 抗 旋转 攻击 鲁 棒 性 平均 提高 了 13.26%， 与 GH 旋 
iE oth Zernike 正 交 线 算法 相 比 ， 抗 旋转 攻击 哩 棒 性 分 别提 高 了 1.1% 和 0.9490; 对 常规 攻击 、 缩 放 攻 击 和 循环 平移 
以 及 小 范围 的 剪 切 攻击 均 具有 较 强 章 棒 性 
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M P) Anti-geometric rotation attack zero watermarking algorithm 
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Abstract: To solve the problems of the zero-watermarking algorithms with weak against geometric attack, this paper proposed 
a zero-watermarking algorithm against geometric rotation attack. Firstly, the algorithm determined the near-destructive 


maximum inscribed square region with approximately pixels lossless in the center area, according to the pixel distortion of the 


he image of the scale-invariant feature transform (SIFT) rotation correction. Then, the square region performed two-level 
redundant discrete wavelet transformation and extracted the low-frequency region. What's more, it extracted the largest 
singular value of each block to construct a transition matrix into block processing from the low-frequency area. Next, to obtain 
a characteristic matrix, it compared the value of each element in the transition matrix with its average value. Finally, it used the 
watermark image and constructed the characteristic matrix a zero watermark. The experimental results show that compared 
with the only rotation correction algorithm of SIFT, the robustness against rotation attack is up by an average of 13.2696. 
Compare with the rotation corrections algorithms of the GH rotation moment and pseudo-Zernike orthogonal moment, the 
anti-rotation attack robustness is up by 1.1% and 0.9496 respectively. Also, it has a good effect on common conventional 
attacks, scaling attacks, cyclic translation and small-scale shear attacks. 

Key words: scale-invariant feature transform; in-square square region; redundant discrete wavelet; Arnold; eigen-matrix; 


singular value decomposition (SVD); anti-geometric attack 


统 数字 水 印 ， 零 水 印 则 能 更 好 地 克服 几何 攻击 。 为 了 有 效 提高 
图 像 水 印 算法 抗 几 何 攻击 能 力 ， 基 于 图 像 特征 点 的 第 2 代 图 像 
零 水 印 中 的 思想 ， 主 要 通过 利用 原始 载体 图 像 的 内 部 特征 水 印 技术 得 到 外 界 了 广泛 关注 站。 文献 [3~5] 利 用 图 像 尺度 不 变 
进行 构造 ， 保 证 了 原始 载体 图 像 完整 性 ， 打 破 了 传统 水 印 算 法 特征 变换 (scale-invariant feature transform, SIFT) 匹配 攻击 前 
鲁 棒 性 与 不 可 见 性 之 间 的 了 矛盾。 一直 以 来 ， 由 于 几何 攻击 会 破 后 的 特征 点 作为 模板 来 校正 图 像 的 几何 失真 ， 从 而 达到 同步 水 
坏 载 体 图 像 与 水 印 的 同步 性 ， 使 得 水 印 无 法 正确 提取 。 相 比 传 印信 息 抵抗 旋转 攻击 的 目的 ， 但 均 没 有 考虑 到 旋转 攻击 带 来 的 
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像素 缺失 问题 。 贾 超 等 人 跨 则 提出 了 基于 改进 SIFT 的 抗 几何 攻 
击 的 水 印 算法 ， 将 旋转 不 变 纹 理 特征 信息 融合 到 传统 SIFT 特 
征 向 量 当 中 ， 进 而 提高 了 匹配 进度 。Zhang 等 人 通过 使 用 加 
速 鲁 棒 特征 (speeded-up robust features, SURF) 和 RANSAC 
算法 整体 提高 了 提取 特征 点 的 质量 ， 再 使 用 放射 变换 来 提高 校 
正 的 精度 ， 其 本 质 是 优化 了 特征 点 检测 ， 并 没有 考虑 旋转 攻击 
带 来 的 像素 值 缺 失 导致 提取 的 水 印 缺失 问题 。 

虽然 特征 点 匹配 算法 能 够 有 效 抵抗 旋转 攻击 ， 但 是 随 着 旋 
转角 度 的 增加 ， 算 法 的 鲁 棒 性 也 将 越 来 越 差 。 其 中 ， 最 直接 原 
因 是 旋转 攻击 不 仅 将 元 素 位 置 进行 改变 ， 同 样 也 会 损失 一 定 的 
像素 值 ， 造 成 旋转 校正 后 的 图 像 一 部 分 的 像素 缺失 。 因 此 针对 
该 现象 , 传统 水 印 算法 也 通过 构建 不 变 域 或 不 变 矩 的 方法 构造 。 
Jia 等 人 外 提出 了 基于 局 部 信息 图 像 归 一 化 的 抗 几 何 攻 击 数 字 
水 印 算法 ， 主 要 通过 选择 一 个 旋转 校正 后 不 会 带 来 影响 局 部 的 
GH 区 域 进 行 水 印 的 供 入 与 提取 ,但 随 着 舱 入 水 印 区 域 的 缩小 ， 
算法 的 容量 也 会 减少 。 此 外 ， 朱 丹 丹 等 人 外 将 水 印信 息 矢 入 到 
载体 图 像 低频 的 伪 Zernike 和 珑 幅 值 当 中 ， 提 高 了 水 印 鲁 棒 性 。 
陈 青 等 人 0 运用 Harris 算法 提取 载体 图 像 特征 点 ， 并 构建 局 部 
特征 区 域 ， 通 过 量化 调制 伪 Zernike 矩 幅 值 将 水 印信 息 骨 入 到 
局 部 特征 区 域 ， 利 用 伪 Zernike 和 矩 的 相位 信息 对 受到 旋转 攻 了 
的 图 像 进行 几何 校正 ， 能 够 有 效 地 抵抗 旋转 攻击 。 
文献 [8,10] 可 得 , 通过 寻找 载体 图 像 不 会 因 受 到 旋转 攻击 
而 像素 值 失 真 的 部 分 来 完成 水 印 的 姐 入 ， 可 有 效 地 提高 SIFT 
几何 校正 水 印 算 法 的 鲁 棒 性 。 通 过 多 次 实验 , 若 在 载体 图 像 的 中 
心 点 处 附近 选取 像素 缺损 范围 最 小 且 校 正 后 图 像 的 内 接 正方 形 
区 进行 零 水 印 的 构建 与 水 印 的 提取 , 可 以 减少 像素 的 缺损 现象 。 
此 外 ， 为 了 弥补 因 减 少 载体 图 像 面积 所 带 来 的 嵌入 容量 减少 问 
题 ， 所 以 本 文选 用 元 余 离 散 小 波 变 换 (RDWT) 在 不 减 小 分 解 后 
子 图 大 小 特性 的 前 提 下 提取 图 像 的 稳定 低频 区 域 00， 并 提取 图 
像 块 的 最 大 奇异 值 构造 特征 和 矩阵， 与 加 密 后 的 水 印信 息 共 同 构 
EFKE. 
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1 ”构建 区 域 的 选择 


1.1 尺度 不 变 特征 变换 

尺度 不 变 特征 变换 (SIFT) 是 用 来 检测 图 像 局 部 特征 的 一 种 
算法 ， 在 多 尺度 空间 659 提取 尺度 、 位 置 、 旋 转 不 变量 中 特征 点 
的 过 程 。 因 SIFT 特征 点 算 子 具有 尺度 不 变性 和 可 改变 图 像 旋 
转角 度 或 亮度 等 特性 ， 所 以 SIFT 变换 的 本 质 内 容 是 构建 尺度 
空间 并 搜索 具备 高 鲁 棒 性 的 图 像 突出 的 信息 特征 点 ， 且 这 些 特 
征 点 不 会 因 光 照 强度 改变 、 遮 挡 等 问题 出 现 而 消失 。 即 在 尺度 
空间 检测 极 值 点 ， 然 后 对 极 值 点 进行 筛选 ， 提 取 在 每 个 稳定 特 
征 点 周围 的 图 像 局 部 特性 形成 局 部 描述 符 。 通 过 欧式 距离 实现 
匹配 ， 得 到 匹配 后 特征 点 的 尺度 、 位 置 、 方 向 和 描述 子 ， 完 成 
图 像 的 几何 校正 过 程 。 
1.1.1 尺度 空间 的 构建 

尺度 空间 的 作用 在 于 模拟 图 像 数据 的 多 尺度 特征 。 高 斯 卷 
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积 核 是 实现 尺度 变换 的 唯一 线性 核 。 用 函数 L(x,y,0) 描 述 图 像 
的 尺度 空间 ， 是 图 像 PG, y) 和 尺度 可 变换 的 高 斯 函数 
G(X,y,0) 卷 积 的 结果 。 即 高 斯 函数 G 对 图 像 了 的 模糊 函数 见 式 
(1) (2)。 


L(x, y, o) = G(x,y,0) 8 I(x, y) (1) 
-Q3 y?) 
G(x,y,0)- ze Jo (2) 
2zo 


其 中 : O 表示 卷 积 ，(Xx, y) 表示 空间 坐标 ， 表 示 图 像 的 位 置 ， 
O 表示 尺度 空间 因子 ， 其 值 的 大 小 与 图 像 了 被 平滑 的 程度 成 正 
比 。 当 图 像 被 平滑 得 越 多 ， 则 图 像 呈 现 得 越 模糊 ， 相 应 的 尺度 
越 大 。 


DoG( 高 斯 差分 尺度 空间 ) 是 通过 将 高 斯 金字 塔 中 每 一 组 中 
上 下 相 邻 的 两 层 高 斯 尺度 图 像 进行 相 减 而 得 。 高 斯 差分 函数 见 
式 (3)。 


D(x, y,0) = (G(x, y, ko) - G(x, y, 0)) G I(x, y) 
= L(x, y, ko)- L(x, y,o) 
1.1.2 关键 点 方向 参数 确定 及 特征 点 描述 算 子 确定 
关键 点 方向 参数 的 确定 主要 利用 边缘 强度 M 和 边缘 方向 
0, SR SR. 


(3) 


Me; Q3) 2 LG. 3) -L-1 y € Qs y D - LG y-D)* 

a (L(x, y+1)- L(x, y-1)) (4) 
(L(x+1,y) - Lix-l y) 

其 中 : L 表示 关键 点 所 在 尺度 。 通 过 利用 直方 图 统计 邻 域内 像 

素 的 梯度 方向 ， 特 征 点 直方 图 的 峰值 可 表示 该 特征 向 量 的 局 部 

方向 ， 再 计算 特征 的 主 方向 ， 保 证 描述 符 具备 旋转 不 变性 。 特 

征 向 量 的 生成 过 程 如 图 1 所 示 。 图 中 黑 点 表示 某 特征 关键 点 ， 

每 个 方块 小 格 表示 关键 点 邻近 点 的 一 个 像素 ， 像 素 的 梯度 方向 

为 箭头 所 指 方向 ， 梯 度 的 模 即 为 箭头 的 长 度 。 


Osrr (x, y) = tan 


X] 1 特征 点 邻 域 梯度 信息 与 特征 向 量 
1.1.3 SIFT 特征 点 匹配 校正 

a) SIFT 特征 匹配 是 凭借 特征 点 之 间 的 相似 度量 所 进行 的 
过 程 。 先 通过 计算 两 幅 图 像 SIFT 特征 点 描述 符 的 欧 氏 距离 对 
SIFT 特征 向 量 进行 匹配 , 在 获取 特征 向 量 值 后 ,采用 优先 K-D 
树 进行 优先 搜索 查找 每 个 特征 点 的 近似 最 近邻 特征 点 。 检 测 两 
个 特征 点 ， 若 较 短 的 距离 除 以 次 短 的 距离 不 高 于 某 个 上 限 值 ， 
则 接受 此 对 匹配 点 。 通 过 更 改 这 个 上 限 值 ， 可 控制 特征 匹配 点 
的 数量 。 上 限 值 变 大 时 匹配 点 数目 增多 , 产生 较 多 的 误 匹 配点 ， 
匹配 稳定 性 变 差 ， 反 之 同 理 。 通 过 大 量 实验 得 出 ， 当 比例 阔 值 
在 [0.4,0.6] 时 ， 匹 配 效果 理想 。 

b) SIFT 是 一 种 基于 空间 的 对 图 像 旋转 、 缩 放 和 平移 保持 不 
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变性 的 图 像 局 部 特征 描述 算 子 。 若 在 中 心 区 域 存在 两 个 点 A, 
B， 若 图 像 旋转 C 角度 ， 则 认为 A. B 两 点 均 旋转 相同 角度 ， 
即 只 需 判 定 两 对 特征 点 之 间 连 线 的 夹 角 实现 旋转 校正 其中， 
oo 为 两 图 像 中 心 点 的 连 线 ) 。 提 取 并 匹配 旋转 攻击 前 后 两 幅 图 
像 SIFT 特征 点 ， 利 用 匹配 特征 点 的 位 置信 息 进 行 旋转 校正 。 
假设 原始 载体 图 像 中 的 两 个 特征 点 坐标 分 别 为 (wyi) 和 


(x,»;) ; 经 过 旋转 攻击 后 的 两 个 特征 点 的 坐标 分 别 为 (%,») 和 
(xoy) ， 旋 转角 度 图 如 图 2 所 示 。 
O D 
Bu) —— BG; Nj ) 
图 2 特征 点 旋转 示意 图 
引 使 用 两 对 特征 点 之 间 连 线 的 水 平 夹 角 的 变化 矫正 图 像 ， 
利用 公式 (5) 求 出 每 组 的 旋转 角度 。 
Ag=arctan -arctan s. (5) 
jj XX; 


b) 将 所 有 分 组 求 出 的 旋转 角度 进行 n 等 分 ， 利 用 直方 图 求 
出 各 个 区 间 的 频数 m. 及 相应 的 旋转 角度 p, 找 到 最 大 频数 mmax， 
求 出 频数 相对 较 大 的 旋转 角度 pm， 第 选 在 [pm-1,pm+1] 内 的 旋 
转角 度 ， 得 出 筛选 后 旋转 角度 的 平均 值 即 为 旋转 校正 角度 值 ， 
见 式 (6)(7)。 


num 
P»? (pCi) x mQ)) 
pm-L z (6) 
2 m) 
i-l 
1 N 
m YA 
iem 之 e, (7) 
其 中 : num 表示 频数 m(i) 大 于 0.85Xmmax[pm-1，pm+1] 区 间 内 
的 旋转 角度 ，N 表示 筛选 后 的 旋转 角度 个 数 。 


1.2 ”内 接 正方 形 安全 区 域 的 确定 
1.2.1 旋转 校正 后 像素 缺损 的 规律 

SIFT 载体 图 像 在 经 过 1”~359” 的 旋转 过 程 中 边缘 信息 会 
有 程度 不 一 的 像素 缺损 ， 且 经 过 旋转 校正 后 的 载体 图 像 也 无 法 
弥补 之 前 缺失 的 像素 部 分 。 原 始 载体 图 像 和 受到 不 同 角度 旋转 
攻击 的 效果 如 图 3 所 示 , 仅 利 用 SIFT 算法 校正 后 的 效果 如 图 4 
所 示 ， 不 同 角度 旋转 攻击 后 的 旋转 校正 结果 以 及 缺损 情况 见 表 
1。 


由 图 4 可 得 如 下 启发 : 
a) 无 论 经 过 任何 角度 的 旋转 ， 原 
有 一 部 分 像素 数据 不 会 出 现 缺损 现象 ; 
b) 载 体 图 像 经 过 15* 的 旋转 攻击 效果 图 与 经 过 
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带 总 


台 载体 图 像 的 中 心地 


75° 的 旋 


转 攻击 效果 图 


c) 当 载体 图 像 受 到 45° 


刘 万 军 ， 等 
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所 缺损 的 像素 


四 积 相同 ; 
s 133° 


. 225° I 315° 


图 像 所 损失 的 面积 达到 极 值 且 相同 (NC 值 )。 


(a) 原始 载体 图 像 


(b) 旋转 5° 


(d) 旋转 45° (e) 旋转 75” (f) 旋转 125* 
图 3 原始 图 像 和 受到 不 同 角度 旋转 攻击 后 的 图 像 


时 ， 载 体 


(a) 校正 旋 


Bro 


(b) 校正 旋转 15? 


(c) 校正 旋转 45? 


(gd) 校 正 旋转 75° 


(e) 校 正 旋转 125° 


图 4 仅 利 用 SIFT 旋转 校正 后 的 图 像 
dl 旋转 攻击 后 旋转 校正 结果 及 缺损 情况 

旋转 角度 (*) ”校正 参数 Lena Boat Bridge ”图像 是 否 有 损失 

5 -5 -4.9972  -49863 -5.0245 有 

10 -10 -10.0120 -9.9973 -9.9876 有 

15 -15 -15.0010 -14.9968 -14.9988 有 

20 -20 -19.9738 -19.9743 — -19.9933 有 

30 -30 -30.0118 -30.0378 -29.9832 有 

45 -45 -45.1030 -44.9863 -44.9899 有 

55 -55 -55.0100 -54.9871 -54.9921 有 

75 -75 -75.0914 -74.9789 -74.9863 有 

由 此 可 得 ， 图 像 若 在 经 过 0~360” 各 个 角度 的 攻击 下 ， 所 
损失 的 面积 以 及 具体 的 位 置 将 会 出 现 周 期 性 ， 即 可 得 
(45K)? ~ (45(K+1)),(k=0,1...7) 规 律 ， 呈 45” 为 一 周期 变化 。 
特别 地 ， 当 载体 图 像 受 到 45”、135”、225” 和 315” 的 旋转 
攻击 时 ， 损 失 的 面积 最 大 ， 且 损失 的 位 置 相同 ， 因 此 在 经 过 旋 
转 校正 后 的 45” 图 像 中 找到 内 接 的 正方 形 ， 其 他 角度 同 理 可 
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a. 1.3， 宛 余 离散 小 波 
因此 ， 本 文 在 SIFT 旋转 校正 算法 的 基础 上 ， 另 在 载体 图 由 于 宛 余 离散 小 波 变换 每 进行 下 一 级 分 解 时 所 得 的 每 个 子 
像 的 中 心 处 寻找 某 一 固定 安全 区 域 ， 使 得 无 论 旋转 角度 大 小 ， “图 矩阵 规模 与 原始 图 形 矩 阵 规模 一 样 ， 所 以 利用 其 特性 不 仅 可 


EE 


该 处 所 提取 的 水 印信 息 在 受 旋转 攻击 等 几何 攻击 的 干扰 下 具有 ”以 提取 到 图 像 稳定 的 低频 区 域 ， 而 且 低频 区 域 的 矩阵 规模 会 与 
较 强 的 健壮 性 。 原始 图 像 保 持 一 致 ， 进 而 提高 算法 的 鲁 棒 性 和 容量 。 二 乡 
1.2.2 安全 区 域 的 确定 RDWT 分 解 示意 图 如 图 7 所 示 。 

在 人 类 视觉 角度 考量 , 当 旋 转角 度 为 (45”、135”、225° 
和 315”) 时 ， 所 选 的 区 域 〈 红 框 标注 部 分 ) 刚好 零 损 失 ; 当 m a NS 
角度 超过 A45" 时 ， 所 选择 的 区 域 同样 安全 ， 所 选择 的 区 域 不 会 
受到 旋转 攻击 而 带 来 损失 ， 故 本 文 将 所 选择 的 区 域 叫做 内 接 正 ud aZ D 
方形 安全 区 域 。 选 取 构 建 区 域 原理 以 及 效果 如 图 5 所 示 。 2d vr un | mm 

(a) 一 级 RDWT 分 解 (b) 二 级 RDWT 分 解 


图 7 二 级 RDWT 分 解 示 意图 
设 图 像 矩 阵 A(Qn X m. 图像 A 经 过 一 级 RDWT 分 解 得 到 一 
个 低频 信息 ，L 万 表示 图 像 的 逼近 部 分 ， 其 矩阵 大 小 为 (ma Xm), 


(a) 校正 旋转 45” — (b) 构 建 正方 形 区 域 (ec) 标 记 所 选 的 正方 形 LHi, HLi, HHi 分 别 表示 图 像 的 垂直 、 水 平和 对 角 线 方向 的 高 
图 5 选取 构建 区 域 原 理 图 以 及 效果 频 信 息 , 三 者 表示 图 像 的 细节 部 分 , 其 大 小 同样 分 别 为 (m Xm). 


针对 旋转 攻击 带 来 图 像 失 真 以 及 旋转 校正 后 图 像 像素 值 缺 。 图 7(b) 则 是 对 LL 进行 下 一 级 的 RDWT 分 解 , 同 理 得 到 三 个 高 
失 的 问题， 文献 [10] 利 用 伪 Zernike 矩 的 重 构图 像 形成 近似 圆 ， 。 频 信息 和 一 个 低频 信息 ， 其 矩阵 大 小 保持 不 变 。 其 中 ， 小 波 分 
但 在 像素 选取 以 及 分 块 处 理 相对 麻烦 ， 故 未 文选 取 内 切 正方 形 。 解 级 数 越 高 ， 所 得 深层 低频 子 带 越 集中 包含 了 被 分 解 图 像 的 绝 
作为 安全 区 域 。 对 于 内 接 正方 形 模 长 的 选取 问题 ，_ 方 面 ， 若 ”大 部 分 信息 。 相 比 高 频 信息 易 受 到 噪声 的 攻击 ， 低 频 子 图 抵抗 
该 正方 形 区 域 的 面积 选取 过 大 ， 会 造成 更 多 的 像素 缺损 ， 另 一 ”外 来 影响 因子 能 力 较 好 ， 其 稳定 性 较 强 。 
方面 ， 若 选 定 面积 选取 过 小 ， 则 会 减少 水 印 的 容量 。 因 此 ， 为 与 离散 小 波 变换 (discrete wavelet transform，DWT) 相 比 ， 
保证 所 提取 的 特征 矩阵 图 像 部 分 不 受到 损失 目 尽 可 能 多 的 包含 DWT 变换 后 所 得 到 一 级 子 图 矩阵 规模 为 原始 图 像 矩 阵 的 1/4, 
水 印信 息 ， 本 文 将 在 受到 损失 面积 最 大 的 旋转 校正 后 的 图 像 中 Ti RDWT 变换 后 子 图 矩阵 规模 不 变 , 因此 在 水 印 算法 容量 角度 
确定 内 接 正方 形 的 规模 并 保持 边 长 度 不 变 ， 使 得 即使 在 其 他 旋 。 上 ，RDWT 更 具有 优势 ， 同 样 文献 [11] 也 印证 了 应 用 在 水 印 算 
转角 度 下 ,该 内 接 正方 形 区 域内 的 像素 点 仍 不 发 生变 化 . HE, 法 中 RDWT 比 DWT 具有 更 强 的 鲁 棒 性 。 因 此 选择 RDWT 变 
DETRE, ESERTEPE — MUNDIAL AUS SVD MON REOR UR RERN 
区 域 具有 一 定 的 可 操作 性 。 棒 性 。 
为 了 验证 所 选 的 区 域 在 受到 旋转 攻击 时 ， 安 全 区 域内 是 否 . 
出 现 像素 数据 损失 的 现象 ， 载 体 图 像 经 过 各 种 角度 旋转 攻击 的 。 2 。 零 水 印 的 构建 和 水 印 提取 
效果 (正常 情况 下 内 切 正方 形 的 规模 ) 如 图 6 所 示 。 本 文 算法 主要 针对 旋转 校正 后 的 相关 水 印 算法 进行 优化 ， 
DW 因此 在 旋转 校正 阶段 在 文献 [5] 的 SIFT 方法 校正 基础 上 ， 另 增 
加 了 内 接 正方 形 安全 区 域 确定 理论 。 算 法 具体 分 为 零 水 印 的 构 
建 与 水 印 的 提取 两 个 步 又 。 在 构建 零 水 印 之 前 需要 对 水 印 进行 
Arnold 置 乱 加 密 ， 以 提高 算法 的 安全 性 和 和 鲁 棒 性 。 水 印 图 像 置 
乱 后 的 效果 如 图 8 所 示 。 


B le is 
q E P ra 
L pL zu MR Au 


(a) 原始 水 印 图 像 (b) WSUSHPKEUER 
À 图 8 水 印信 息 置 乱 前 后 的 效果 
(d) 旋转 45" (e) 旋转 55" (D 旋转 75。 2.4. 零 水 印 的 构建 流程 


图 6 旋转 攻击 后 选择 区 域 缺损 情况 a) 对 水 印 图 像 进行 Arnold 置 乱 ， 得 到 加 密 后 的 水 印 图 像 
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W(Q(nxn) 并 保存 秘 钥 K- 

b) 根据 1.2 节 的 讨论 ， 确 定 载体 图 像 中 不 会 因 受 旋转 攻击 
而 导致 像素 缺失 的 内 接 正 方形 MGnxm)， 其 中 mod(m,n)-0, 
保存 内 接 正方 形 的 左上 角 坐 标 以 及 边 长 。 
c) 对 M 进行 二 级 RDWT 变换 ， 提 取 其 逼近 子 图 并 分 块 
B, P, ij=1,2,...m/n. 

d) 对 图 像 块 进行 SVD 分解 , 并 提取 每 个 块 的 最 大 奇异 值 ， 
构建 过 渡 和 矩阵 了 。 

e) 根据 了 中 元 素 值 与 其 平均 值 mean_7Y 的 大 小 关系 根据 特 
征 矩 阵 C， 如 式 (8) 所 示 。 


1 Y(j)2mean Y 
C(,j)- A (8) 
0 Y(i,j)<mean_Y 


PRERE C 与 加 密 水 印 W HEITER, ESRUS 
息 零 水 印 CW， 并 送 至 第 三 方 权威 机 构 认 证 中 心 IPR) 注 册 。 
2.2 水 印 的 提取 

a 人 对 可 能 受到 旋转 攻击 的 
待 检测 图 像 P。 

b) 按 照 构建 零 水 印 的 步骤 b)~e， 同 理 得 到 特征 矩阵 C*。 

co) 从 第 三 方 权威 机 构 认 证 中 心 获取 零 水 印 CW,— E REGERE 
阵 进行 异 或 操作 ， 并 提取 秘 钥 K 进行 逆 Arnold 置 乱 ， 得 到 提 
取 的 水 印 图 像 W*。 


R] 


像 进行 旋转 判断 与 校正 ， 得 到 


真实 验 采 用 MATLAB 2014a 作为 实验 平台 , 选取 6 幅 纹 
理 不 同 的 灰 度 图 像 (512X512)， 水印 图 像 为 32X32 的 带 有 “过 
宁 工 大 ”标志 的 二 值 图 像 。 其 中 置 乱 秘 钥 K 为 20。 本 文采 用 误 
人 码 率 (bit error ratio, BER) 和 归 一 化 相关 系数 (normalized 
cross-correlation, NC) 来 衡量 提取 到 的 水 印 图 像 与 原始 水 印 
像 之 间 相 近 程 度 。BER 是 指 提取 的 水 印 错误 的 比特 数 占 水 印 全 
部 比特 数 的 比例 , 这 个 值 介 于 0~1 之 间 , 显然 BER 越 接 近 于 0， 
水 印 的 鲁 棒 性 越 好 。NC 值 的 范围 在 0-1 之 间 ， 且 该 值 越 接 近 
于 1， 表明 该 数字 水 印 算法 的 鲁 棒 性 越 好 。 

BER 定义 为 

BER =N, / N,,, x 100(96) (9) 

其 中 : NMemor 表 示 与 原始 水 印 相 比较 ， 提 取 的 水 印信 息 中 错误 的 
ERMER AZG Nois 表示 水 印 中 总 的 比特 数 。 

NC 定义 为 


» (10) 
2,6 j} 
其 中 ， wi 四 为 原始 水 印 图 像 ，w(i,j) 为 提取 的 水 印 图 像 ，M 
为 水 印 图 像 的 长 或 宽 。 
3.1 旋转 攻击 测试 

本 节 仅 从 以 原始 载体 图 像 的 中 心 点 为 旋转 中 心 进行 旋转 攻 
击 测试 。 对 受到 旋转 攻击 的 载体 图 像 ， 先 进行 SIFT 特 


征 点 校 J 


所 示 ， 从 不 同 


6 ilii [e] 


图 像 抵抗 
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区 


载体 图 像 受到 不 同族 
图 像 仍然 清晰 可 见 ， 而 且 从 肉眼 训 
旋转 攻击 的 NC 值 数 扩 
体 可 得 ， 当 受到 任意 角度 旋转 攻 了 
均 在 0.996 2 以 上 ,尤其 当 旋转 90” 时 ,其 NC 值 均 为 1。 原 
在 于 当 图 像 受 到 90° 旋转 攻击 时 ， 图 像 无 任何 像素 值 缺 损 ， 
有 旋转 不 变性 ， 所 以 90° 的 旋转 攻击 对 水 印 的 提 : 
x ses 180° 30270? 的 旋转 攻击 


H. SVD 
无 影响 ; 同 理 ， 当 
取 到 的 NC 值 也 为 1。 


FE 方法 ， 然 后 选 定安 全 区 域 ， 从 而 提 
测试 实验 将 对 6 个 不 同 的 

测试 图 像 受到 旋转 攻击 效果 和 提取 的 水 印 
Jt Pr ff RE Bid 


it 可 以 进行 


Fi, 


其 提 : 


取 


Cape EAR. 
出 水 印信 息 。 本 
像 进行 旋转 攻击 测试 。 
图 像 如 图 9 
UI 
版 权 的 认证 


区 的 水 印 


o 


EILE 2。 从 表 2 整 
取 到 的 水 印 NC 


辽宁 
工大 
(e) Bridge 旋转 45? (f) Couple 旋转 45° 
图 9 测试 图 像 受到 旋转 攻击 效果 和 提取 的 水 印 图 像 
表 2 不 同 图 像 抵 抗 不 同 角 度 旋转 攻击 的 NC 值 ( 逆 时 针 旋 转 ) 
旋转 角度 () Lena Barbara Baboon Boat Bridge Couple 
1 1.0000 0.9993 0.9987 1.0000 1.0000 0.9993 
5 0.9974 0.9980 0.9993 0.9980 0.9993 0.9980 
10 0.9974 0.9975 0.9955 0.9980 1.0000 0.9974 
15 0.9980 0.9968 0.9955 0.9955 0.9993 0.9980 
25 0.9993 0.9981 0.9948 0.9961 0.9987 0.9974 
35 0.9980 0.9962 0.9929 0.9955 0.9993 0.9974 
45 0.9974 0.9974 0.9942 0.9974 0.9993 0.9968 
55 0.9993 . 0.9968 0.9968 | 0.9961 0.9993 0.9974 
65 0.9980 | 0.9962 0.9942 0.9987 0.9980 0.9968 
75 0.9980 | 0.9962 0.9968 0.9974 0.9993 0.9968 
85 0.9980 0.9987 0.9974 0.9993 1.0000 0.9993 
90 1.0000 1.0000 1.0000 1.0000 1.0000 1.0000 
120 0.9980 0.9987 0.9935 0.9980 0.9987 0.9968 
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3.8 对比 实验 测试 
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常规 的 攻击 主要 有 几何 攻击 与 非 几何 攻击 ， 而 在 3.1 节 中 为 了 进一步 体现 本 文 算 法 抗 旋转 攻击 的 优越 性 ， 与 文献 
已 经 讨论 过 本 文 算法 抗 旋转 攻击 的 能 力 ， 因 此 在 本 节 中 不 进行 [5,8,13] 进 行 对 比 。 其 中 文献 [3] 只 采用 了 SIFT 进行 校正 ， 却 没 
旋转 攻击 测试 。 取 6 幅 常规 图 像 中 的 Lena, Barbara, Baboon 有 考虑 旋转 攻击 给 旋转 校正 后 图 像 带 来 的 损失 ， 而 后 两 种 算法 
为 例 ， 其 抵抗 非 几何 攻击 测试 结果 见 表 3。 也 仅 考 虑 到 了 旋转 却 没 有 校正 的 问题 ,文献 [8] 采 用 
从 表 3 中 的 数据 可 得 ， 随 着 攻击 强度 的 增加 ， 提 取水 印 的 ” Gaussian-Hermite 确定 局 部 不 变 区 域 构建 特征 矩阵 ， 文 献 [13] 
NC 值 逐 渐 减 小 , 但 减 小 速度 缓慢 , 其 整体 最 小 NC 值 也 在 0.95 ”通过 计算 图 像 归 一 化 后 的 伪 Zernike 和 矩 ， 然 后 选取 部 分 合适 的 
以 上 。 特 别 是 针对 JPEG 压缩 攻击 ， 当 JPEG 压缩 的 质量 因子 ”和 拖 通 过 量化 调制 嵌入 水 印信 息 。 从 水 印 算法 方面 ， 文 献 [5] 选 择 
为 1 (压缩 比 为 100:1000 ， 从 三 幅 不 同 的 图 像 提 取 到 的 水 印 ”掩蔽 性 好 的 图 像 块 作为 嵌入 子 块 ， 得 到 子 块 二 级 离散 小 波 低频 
NC 分 别 为 0.971 6, 0.972 3 和 0.951 0, 可 见 本 文 算法 具有 抵抗 的 左 奇异 和 矩阵， 通过 修改 第 一 列 元 素 间 的 大 小 关系 来 实现 水 印 
高 强度 JPEG 压缩 攻击 的 性 能 。 同 理 ， 本 文 算法 对 其 他 类 型 攻 WRA: 文献 [13] 计 算 图 像 归 一 化 后 的 伪 Zernike 3B, 然后 选取 
击 的 效果 也 较为 显著 ， 均 近似 可 得 1。 部 分 合适 的 和 矩 通过 量化 调制 谋 入 水 印信 息 ; 而 本 文采 用 非 嵌 入 
幅 图 像 受 到 不 同 几 何 攻击 后 提取 水 印 的 数据 见 表 4。 其” ” 式 的 零 水 印 方法 ， 不 会 给 载体 图 像 带 来 任何 的 视觉 影响 ， 因 此 
中 ， 由 于 采用 SIFT 特征 点 的 几何 校正 ， 使 得 算法 能 够 有 效 地 ”本 文具 有 较 强 的 不 可 见 性 。 不 同文 献 方法 下 抗 旋转 攻击 的 对 比 
抵抗 缩放 攻击 和 行列 偏 移 攻击 ， 而 且 提 取水 印 NC 值 均 接近 1, 实验 数据 见 表 5. 
T 具备 良好 的 抗 几何 攻击 能 力 。 当 受到 剪 切 攻击 时 ， 由 于 剪 切 载 表 5 不 同文 献 的 旋转 攻击 对 比 实验 测试 结 
E 体 图 像 的 1/64 时 , 没有 对 提取 图 像 特征 点 的 区 域 造成 影响 ， 进 旋转 攻击 角度 /" ”文献 [5] — 文献 [8] — 文献 [13] 。” ”本文 
v ”而 不 会 影响 特征 点 的 提取 ， 从 而 使 提取 的 水 印 非常 完整 ， 而 当 0 1.0000 1.0000 1.0000 1.0000 
Lam 增加 剪 切面 积 时 ， 会 把 部 分 提取 特征 点 的 区 域 减 掉 ， 直 接 影响 10 0.9366 0.9867 0.9930 ^ 0,9974 
É 了 提取 水 印 的 质量 ,因此 本 文 算法 能 够 抵抗 小 范围 的 剪 切 攻击 。 20 0.8966 0.9820 09710 0.9993 
O 表 3 抵抗 非 几 何 攻击 测试 结果 (CNC 值 ) 30 0.8446 0.9855 0.9830 0.9980 
© 攻击 方式 参数 Lena Barbara Baboon 40 0.7325 0.9822 0.9860 0.9987 
- 高 斯 噪声 Oue 02211 02303 02863 50 0.7339 0.9820 0.9980 0.9980 
e i ud d dn 60 0.8447 0.9873 0.9890 0.9961 
4 a 椒盐 噪声 Mio boh Mad 70 0.8945 0.9865 0.9800 0.9987 
~ 0.020 0.9923 0.9929 0.9923 
,一 i 09716 09723 0.9510 80 0.9087 0.9898 0.9880 0.9947 
X JPEG 10 0.9974 0.9929 0.9832. 90 1.0000 0.9901 1.0000 1.0000 
E 20 0.9936 0.9974 0.9916 由 表 5 中 的 数据 可 分 析 ， 通 过 旋转 校正 优化 后 的 水 印 算法 
C etx MESE Re 抵抗 旋转 攻击 的 鲁 棒 性 要 明显 优越 于 只 进行 校正 的 水 印 算法 ， 
O BE Ne 其 中 文献 [8,13] 和 本 文 算法 属于 旋转 校正 后 选择 合适 的 区 域 进 
mmama SO CHER 0990 RBS 行 水印 的 提取 ,而 文献 [] 只 进行 旋转 校正 。 文献 [8] 抵 抗 旋转 攻 
x5 0.9955 0.9968 0.9955 | M m 
"€ 3x3 0.9961 0.9987 0.9980 A NC E Bi ed plo NC 
5x5 0.9948 0.9955 0.9942 值 为 0.988 8， 本 文 抵抗 旋转 攻击 的 平均 NC (873 0.998 1, 可见 
图 像 锐 化 一 0.9916 0.9833 0.9871 本 文 算 法 抵抗 旋转 攻击 的 性 能 要 略 强 与 文献 [8,13]。 本 文 算法 与 
直方 图 均衡 化 Ex por 0:3872 0.2838 SCEA[S]2S UL, ABAE RE R K RKA A XCTI AEK ER EE 
表 4 不 同 图 像 抵抗 几何 攻击 测试 结果 (NC 值 ) 本 文 算法 优 于 文献 [8] 的 主要 原因 就 是 采用 了 图 像 矩 阵 的 奇异 
攻击 方式 参数 Lene Barbara Baboon 值 ， 图 像 矩 阵 的 奇异 值 能 够 抵抗 扰动 干扰 ， 能 够 有 效 抵 抗 旋转 
Hes Dies HERE 02243 校正 后 ， 通 过 插值 方式 进行 图 像 复原 的 干扰 。 而 文献 [13] 虽 然 
缩放 攻击 03 02280 0.3287. 0:22 能 够 在 一 定 程度 上 改善 抵抗 旋转 攻击 的 鲁 棒 性 ， 但 其 水 印信 息 
Hm 1:0009 1:0000 10000 仅 为 长 度 为 128 的 无 任何 意义 的 二 值 序列 ， 且 水 印 算法 容量 较 
ma 左上 角 1/64 1.0000 1.0000 1.0000 小 。 
剪 切 攻击 
左上 角 1I6 09871 09878  — 0.9656 KO 常规 攻击 对 比 实验 提取 水 印 的 NC fi 
(50,50) 1.0000 1.0000 1.0000 攻击 方式 参数 文献 [5] 文献 [13] 本 文 算法 
— (100,0) 1.0000 1.0000 1.0000 n 9b (sn — mem 
(0,100) 1.0000 1.0000 1.0000 is € P T 
(128, 128) 1.0000 1.0000 1.0000 
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中 值 滤波 3x3 0.9937 1.0000 0.9968 
攻击 5x5 0.9894 0.9670 0.9929 
0.01 0.9464 0.9500 0.9911 

高 斯 噪声 
0.02 0.8522 0.8400 0.9852 
0.01 0.9623 1.0000 0.9968 

椒盐 噪声 
0.02 0.9350 0.9490 0.9923 
0.5 0.9988 0.9693 0.9981 

缩放 攻击 
1.5 1.0000 1.0000 1.0000 
(50,50) 1.0000 1.0000 1.0000 
(100,0) 1.0000 1.0000 1.0000 

循环 平移 
(0,100) 1.0000 1.0000 1.0000 
(128,128) 1.0000 1.0000 1.0000 
剪 切 攻击 左上 角 1/4 0.8661 0.9783 0.8488 
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: 一 种 抗 几何 旋转 攻击 零 水 


本 文 算法 在 空域 上 进行 嵌入 区 域 的 选择 , 当 受到 剪 切 攻击 时 ， 
有 可 能 把 所 选择 的 区 域 减 掉 ， 使 得 提取 水 印 的 区 域 缺失 ， 造 
提取 水 印 的 不 完整 ， 但 当 受 小 面积 的 剪 切 攻击 时 ， 本 文 算法 
JE fe de PLUS ELE RO E P 


Ix] 


HT 


o 


结束 语 


针对 几何 旋转 攻击 导致 图 像 部 分 像素 信息 缺损 无 法 提取 完 
水 印 的 问题 ， 本 文 提 出 了 一 种 抗 几 何 旋 转 校正 后 的 图 像 优 化 
水 印 算法 。 原 始 载体 图 像 预 处 理 阶段 ， 在 SIFT 法 几何 校正 
， 通 过 构建 内 接 正方 形 安全 区 域 ， 实 现 了 在 保证 特征 矩阵 中 
含 像素 信息 无 缺损 下 完成 零 水 印 的 构建 ， 恢 复 了 水 印 的 
。 在 零 水 印 的 构建 和 水 印 提取 阶段 , 利用 RDWT 变换 保留 与 
始 图 像 矩阵 规模 同样 大 小 的 子 图 特性 ， 提 高 了 水 印 算法 的 容 


可 
str 


t 


fi 


; 并 利用 低频 信息 和 SVD 的 稳定 特性 来 提高 算法 的 鲁 棒 性 。 


对 载体 图 像 进 行 JPEG 
以 Lena 图 像 为 例 


10 所 示 。 


0.95 


压缩 攻击 、 噪 声 攻 击 和 滤波 攻击 。 
， 本 文 算法 与 文献 [3,13] 对 受到 不 同 攻击 后 提 
取水 印 的 NC 值 见 表 6。 不 同文 献 旋转 攻击 对 比 实 验 结果 如 图 
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10 不 同文 献 旋转 攻 


表 6 中 可 得 , 本 文 算法 能 


同 程度 JPEG 


Ff 对 比 实验 结果 

了 效 地 抵抗 JPEG JE Hä d 
中 值 滤波 攻击 ， 高 斯 噪声 攻击 、 椒 盐 噪 声 攻 击 、 
环 平移 攻击 和 剪 切 攻击 。 从 受到 不 


Ht 


` 


缩放 攻击 、 循 
压缩 攻击 后 提 


取水 印 情 况 来 看 ， 本 文 算法 提取 水 印 的 NC 值 的 最 小 值 仅 为 
0.985 8， 最 大 值 为 0.998 7， 而 文献 [5,13] 提 取水 印 的 NC 值 的 


最 小 值 分 别 分 0.968 5 和 0.907 0， 
算法 
着 加 
加 突出 。 尤 


lim 


KPE JPEG 压缩 的 能 力 衰减 严 寻 
E JPEG 压缩 攻击 ， 本 文 算 法 抗 JPEG 压缩 攻击 的 能 力 更 
其 是 抵抗 噪声 攻击 , 当 受 到 0.02 的 高 斯 
本 文 算法 的 鲁 棒 性 比 文献 [5,13] 分 别 高 


当 受 到 0.04 的 椒盐 噪声 攻击 时 ,本 


三 


Hx 


hu 


大 值 均 为 1， 


表明 文献 [13] 


， 其 次 是 文 


献 [5]， 可 见 随 


噪声 攻击 时 ， 


H f 15.61% 和 17.2996; 


文 算法 的 鲁 棒 1 


生 比 文献 [5,13] 


分 别 高 出 了 9.95% 和 39.10%。 同样 从 受到 中 值 滤 波 攻击 的 数据 


来 看 ， 本 文 算 法 抵抗 中 值 滤波 攻 了 


f 的 能 力 


也 要 优 于 文献 [5,13]。 


对 于 缩放 攻击 与 循环 平移 攻击 ， 本 文 算法 性 能 与 文献 [3,13] 类 


似 。 但 抗 剪 切 攻击 能 力 ， 本 文 算法 不 如 文献 [5,13]， 其 主要 原 医 


通 
面 


= 


过 大 量 实验 表明 , 算法 能 够 有 效 地 抵抗 旋转 攻击 、 HEROS 
环 平移 攻击 以 及 剪 切 攻击 ， 均 体现 了 良好 的 鲁 棒 性 。 但 在 大 
积 地 实现 剪 切 攻击 方面 ， 本 文 算法 能 力 仍 有 待 提高 。 
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